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IF FUSION BEHAVIOR OF SUPERSATURATED BINARY SOLUTIONS 


Jacob Mazur and Turner Allfrey, Jr. 
Polymer Research Laboratory, The Dow Chemical Company, Midland, Mich. 


PART I 


The nonideality of binary mixtures may be attributed to a number of 
tors, including differences in size and shape. We restrict our discus- 
m to the nonideality arising entirely from the unequal intermolecular 
ces between the molecules. We assume a lattice filled with a mixture 
2 species of spherical molecules, A and B, having identical diameters 
different chemical constitutions. The coordination number of the 
tice will be designated as z. The deviation of the binary mixture from 
tality is determined by the average potential energy of an A-B pair, 


€aA + €BB 


2 


Ae =€ap - 


ere «4, is the energy of association of an A-A pair, egp is the energy 
association of a B-B pair, and e4, is the energy of association of an 
B pair (44> €p8, and egg are negative in sign). Negative values of 
mean that the mixing is exothermic (strong affinity of A molecules for 
nolecules). Positive values of Ae mean that the mixing is endothermic. 
\e is sufficiently positive, limited miscibility results. 
[he statistical thermodynamics of this model have been analyzed by 
ny investigators.* Even the crudest treatment (usually designated as 
‘eroth-order’’ approximation) leads to a reasonable qualitative descrip- 
n of the behavior of such mixtures in spite of a nonzero Ae and the 
sumption of random mixing. 
[he free energy versus composition curve (with the assumption of 
dom mixing) has, for sufficiently large positive values of Ac, 2 minima 
i a central maximum (FicuRE 1). Any homogeneous solution having a 
nposition between n, and n, (the 2 minima) will be unstable relative to 
aixture of 2 separate phases having the compositions n, and n,.If a 
nogeneous solution of intermediate composition n (n, <n< nj) is 
ned, the free-energy curve requires that it should separate into 2 such 
ises. However, there are 2 different kinds of unstable solutions: 
fastable solutions and absolutely unstable solutions. In metastable 
utions (0?F/dn?)r,p > 0 (n is the mole fraction of molecules 4A). 
erefore, minor fluctuations in concentration tend to increase the free 
tgy of the system. The fluctuation in concentration tends to be wiped 
*For an up-to-date review of the problem of the stability of binary mixtures, see 
rogine. 1 
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FIGURE 1, Configurational free energy of binary mixtures. 


out, and the solution, in the absence of seeding, may persist for a 
long time.* In the regions where (07F/dn”)y,p < 0,the F versus n cu 
is concave downward. Hence, any minor fluctuation in concentration le 
to a decrease in free energy; therefore, such a solution is absolu 
unstable with regard to concentration fluctuations, and it will split s 
taneously into 2 different phases. 

The behavior of a solution in the shealutely unstable region. can | 
described. in terms of a ‘‘negative diffusion coefficient,” since a con 
tration gradient will.cause diffusion up the gradient, rather than. dow 1 
and thus accentuate the difference in concentration. . __ 

We propose to reexamine the lattice model in. the statistical. ‘then 
dynamic treatment. described above, to. introduce an elementary diffus 
process (consisting of molecular interchanges), and to. deduce. -from_ 


expanded model an approximate expression for the diffusion coeffici 
as a function of concentration. 


*The possibility of a phase separation in a metastable solution in the absenc 
seed (the so-called heterophase fluctuation) is discussed by Frenkel,2 
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Toward this end we assume the existence of a 1-dimensional concen- 
ation gradient in the positive x direction. In any layer of the lattice 
ented normal to the x direction we assume a constant composition and 
tandom arrangement of the 2 types of molecules. Furthermore, we as- 
ume that n(x) is a continuous function with continuous derivatives. (The 
ase with discontinuities in n(x), that is, when the equilibrium configura- 
ion involves sharp phase boundaries, is discussed below.) 

Let us now consider a particular molecular pair interchange in a binary 
attice. An A-B pair in arrangement will be ‘‘rotated’’ to the interchanged 
trangement II with a certain rate that will depend upon the activation 
nergy E,,;; for the process. The number of A-B interchanges of this 
atticular variety occurring per second across a given plane of the lattice 
ill be proportional to: (The number of sites per layer) - (the probability 
lat a given pair of sites accommodates an A-B pair) - (the probability 
tat the ‘‘atmosphere’’ of neighbors is the particular one) - (eS st Mie y . 
y¥rou/*T), The number of transitions of the precisely reverse form 
II) will be given by a similar expression: (sites/layer) - (probability 
lat a given pair of sites accommodates a B-A pair) - (the probability that 
le ‘‘atmosphere’’ of neighbors is the particular one) - (As /ky 
EtioI/kT | A St is the entropy of activation of the molecular interchange. 
The net flux of molecules A to the right, across the boundary, resulting 
om the interchanges of these 2 detailed types is equal to the difference 
Stween the 2 terms above. It is assumed that AS‘ is the same for all 
etailed molecular interchanges. Although the absolute values of E;5;; 
ad E,;,,; are indeterminate, the difference between them is fixed by m, 
ie net change in the number of A-B ‘‘bonds’’ between the 2 states 
sIGURE 2). 

We now choose a new reference point in such a way that 


> 


) 
; eEruV/kT = gE /kT , g-mAc/2KT 


‘ 
er Erer)/kT _ gE /kT , emNe/2kT 


hus, since we are interested only in the net flux up or down the concen- 
ion gradient, we may factor out not only the indeterminate value of 
4S#/* but also the indeterminate part of each activation energy e = /*7, 
s well as the number of sites per layer, into 1 comprehensive indeter- 
inate constant, K. 
In the absence of a concentration gradient, the above pair of specific 
ementary interchanges will lead to a net flux (because of the asymmetry 
the atmosphere), but this will be just counterbalanced by a net flux in 
e opposite direction by interchanges of a similar type with reversed 


4 
3 


256 Annals New York Academy of Sciences 


POTENTIAL ENERGY 
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FIGURE 2, Change in energy for a typical molecular interchange. 


atmospheres. This offsetting will lead to an over-all net flux of zero. 1 
the presence of a concentration gradient, the over-all net flux is not z 
In some concentration ranges the flux per unit gradient is positive 
cogradient); in some it is negative (or contragradient). The dependen 
the net flux on the concentration will be determined next. 


PART II 


Let us consider an ‘‘atmosphere’’ of N4 molecules A, N-N, molecule 
B, which are distributed among WN sites. Let these sites be located alo 
4 equivalent lines perpendicular to the x axis, as in FIGURE 3. 

Let the coordinates of the neighbors be 0, Ax, 2 Ax, and 3 Ax. 
again assume that the molecules A and B are randomly arranged along th 
lines perpendicular to the x axis. Therefore, the probability that a m 
cule A will occupy 1 of the sites along the line x = 0 is simply Na, 
for brevity, n. For the molecule B, this probability is 1— n. Owing to 1 
concentration gradient along the x axis, this probability will be alt 
as we move along it. Thus, along the line with x = Ax, this probability 
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FIGURE 3. Location of N sites (see text), 


For molecules A: 


n 
n+Ax— 
“~~ * 
For molecules B: 
oe a 
3 l-n-Ax aoe 
Ox 


Let there be a, ‘‘atmosphere”’ sites along the line x =0,a, sites along 


: line x = Ax, and other such sites. Obviously, 3 a =N. Let there be 


molecules A along x = 0, j, molecules A along x = Ax, and soon. 
en the probability C of having an atmosphere of specific structure is: 


om - &=3 j ae 
4 ly, oe iL - ip 
D = / / Lape £2 1-n—-fLAx — 

e Ox re] 


bn, ag ,1) 
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In our diagram a) =a, =1, a, =a,=2; fo=i2=Gf=i3=h 
This “‘probability product’? C can be simplified if all powers of Ax - 
higher than the first are neglected. The simplification is justified, si 


Ax ga << n, everywhere. Then, 


= On = =0 
Can 4 (iaay 4 < de de eS 


It is necessary to multiply this term by a probability that the central pe 
of sites accommodates an AB pair. Since the AB pair is located in 
geometrical center of the atmosphere on the lines £ = 1 and Le- 23 
probability, denoted as D,4z, is given by 


Dap = (n+ ax) (1-n-2hx =) du. 


Multiplying C and D,4p and neglecting again the second powers of axe 
we obtain: 
(Probability of an AB pair) - (probability of the particular 
atmosphere described by Jo > eeeeeesg js) = 


3 
ANA (l= n)V-%4 n(l—-n)+Ax = [a-9 > hi, - 
£=0 


Similarly, for the reverse configuration we obtain 


CDp, =n'4 (l=n)V-NA {x (lee n)+Ax 


x 


P| 

I, 
ay bmy~i) -@n-0]} 4 
2=0 | 


[a-9 Eh-a ys day i) -(on-2)]} 
Z=0 h=0 ! 
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The net flux of A to the right for this particular pair of reverse inter- 
anges is proportional to CD4n eV *7 _ CD, 46 Ue! Tether 
e, the net flux for this interchange is 


@=Kn'4 (l—a)""4 \: (Ton) + Ax = 


n 

. ‘ite ack Solas a as" ou 

[a-a) > hip — 1 oS I (ap i,) -@n-1)] 2" (1,6) 
=o b=o 


eEpsyy/kT 9 -Eyp9q/kT dn -g£ 
ie Ill ig tre’ | = gure kt 
x 


ere K is an indeterminate constant. 


mAe 


Since fag = » Eq» ,=E' —- = , and since E’ is 


Sumed to be independent of the structure of the atmosphere, the net 
x of A to the right for the particular pair of interchanges is 


® =Kn™4 (1-n)"*4 ts (l-n)+Ax = 
x 


3 3 
[a -2) > £j, — nS (ag— jp) - Bn - | a,7) 
g=0 L=0 
sa 2 nial ese — Ax 2 ente/aar} 
J 2kT ox 


As the next step, we shall find the connection between m, the net 
ange in number of AB ‘“‘bonds,’’ and the arrangement of the molecules 
the particular atmosphere. The quantity m is independent of the way 
2 molecules distribute between the sites with L = 0 and $= 1, and sim 
tly, between the sites with y= 2 and aes 3. However, m depends upon 
> way the molecules distribute between the 2 halves of the atmosphere. 
Let there be i molecules A at the left side of the atmosphere, and 
_-i at the right side of the atmosphere, with i = jg+j1; Na-i=jatis: 
erefore, m = 2 (Na — 2i). 

The total flux resulting from all possible interchanges AB and BA is 
sum of EQuATION II,7 over all possible arrangements of the mole- 
The derivation of the total flux from zquaTion II,7 is shown in 
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APPENDIX I. The equation for the total flux is solved explicitly for f 
following cases: 

Casel: z=4; ag =a3=1, a, =a, =2 

Casell: z=6; ag =a, =1, ay =a, =4 

Case Ill: z=8; ag =a; =1, a; =a, =6 


The various terms in the equations for the total flux are arranged in 1 
following way: 


Po ite f2a5 ( Ae: 
Total flux = K [= A, (n) sinh rr -— B, (2) cosh oT 


where A, (n) and B, (n) represent polynomials in n of the order 2 (z—- 
Since A, (n) and B, (n) are always positive, the first term of the ri 
hand side of EquaTion II,8 represents the contribution of the ‘‘co 
gradient’’ flow, while the last term represents the contribution of | 
‘‘cogradient’’ flow toward the total flux, when Ac > 0. When Ac < 0,6 
terms will be cogradient, and we have a case in which the net free ene 
of activation will ‘‘reinforcé’’ the diffusion flow from higher to low 
concentrations. 

It will be worthwhile to show the solution of EquaTion II,8 wh 
Ae¢ =0, as a check of our derivation. In this case, 


ze-l Ze 


1 oN 
Total flux = —K’ >a B, (n) =-K' > 2 a n)\ "NA 


N/2 (N° N 
Koay (7)#a-08 = -K' 
t 


t=0 


ple Sta gs my i — pores 
+ > — 7 


by binomial expansion. 
A further comment should be made about the indeterminate constant J 
This constant includes not only eE’/ kT but the concentration gradie: 


Ax 2 as well. In other words, the net flux is taken now as per ui 


gradient. 


The results of numerical computations are given in FicuREs 4, 5, 2 
6, in which the net flux (expressed as —D/D,) is plotted as a function 
n for several values of A¢/kT. We notice that there exists a critic 
value of Ae below which the flux is always cogradient. This value 
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MOLAR FRACTION n 


FIGURE 4, The dependence of the net flux on concentration for z = 4. 
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FIGURE 5. The dependence of the net flux on concentration for z =6, 
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10 


0 0.2 0.4 0.6 0.8 1.0 
MOLAR FRACTION n 


FIGURE 6, The dependence of the net flux on concentration for z = aa 


Ac (= Aé*) depends not only upon the coordination number z but « 
upon the way the sites are arranged along the diffusion coordinate ( 
other words upon a). The dependence of Ac* (and of D/D,) on the geo 
etry of the lattice arrangements is in striking contrast to the implication 
of the thermodynamic model, where D/D, and Ae* depend only on z : i 
given value of Ac/kT. 


iN adit 


- Part III 


The thermodynamic theory of regular solutions enables us to consti! 
the excess free-energy curve as a function of concentration for conce! 
trated binary solutions. As was pointed out in PART I, these free-energ 
curves in some cases possess 2 ‘‘inflection points.’’ The concentratios 
enclosed by these inflection points (symmetrically arranged aroui 
n = 0.5) correspond to an absolutely unstable region that can be describe 
by negative diffusion coefficients (for example, the parts of FIGURES ‘ 


5, and 6, where the total net flux is greater than zero). These a 
points are defined by 


dn? 


(at given T and p) 
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vese points should correspond to the zero net flux. We shall therefore 
mpare our conditions for zero flux with the solutions of EQUATION 
,1. This is done as follows. 

The statistical thermodynamic theory of regular solutions, based upon 
“zeroth-order’’ approximation, gives the following equation for the 
cess free energy of binary solution (in terms of mole fraction n); 


Ww 
=| .F' 
i 

(III, 2) 


WX and 3 are the number of molecules of the species A and B. F' in- 
udes all terms that enter into the free energy, but that are independent 
| n. The quantity wap is one half of the total increase of potential 
ergy of molecular interchange. From equations III,1 and II,2 we 
ave as a condition for zero net flux 


e F 
ie4Ne - eed [# toga +(-m) log (1 — n) +n (1—n) 


WaB 1 
ae ee a Ill, 
kT 2n (l—n) holt 
? 
Ae 1 
<= ——_— (III,4) 


kT 2zn (1 —n) 


Ficure 7 shows the solution of EquaTion III,4 for z = 4, 6, and 8, 
iere n is given asa function of Ae. In the same plot the conditions for 
so flux, as obtained from ricurREs 4, 5, and 6, are plotted in similar 
shion. We notice good qualitative agreement between the results based 
| our kinetic model and the thermodynamic model. We also notice that 
e conditions for zero flux based on the kinetic model show a stronger 
pendence on z than the thermodynamic model. That is, in the kinetic 
adel increasing z will shift the values of n for which the total flux is 
ro toward terminal points more rapidly than in the case of the thermo- 
lamic model. This difference is not surprising since, as stated before, 
uation for flux based on the kinetic model is a polynomial inn of the 
der of 2 (z — 1), instead of the second order in n, as in the case of the 
ermodynamic model. Also, in our particular models we assumed that 
= 1. Should we reshuffle the neighbors so that ag no longer is unity, 
: obtain a different flux-versus-concentration curve. Therefore, we would 
tain different results for the conditions of zero flux. 
We next derive an eeeretion for the net flux using an expression ob- 
ined by Stearn et al.* and the thermodynamic equation for excess free 
ergy of randomly mixed concentrated binary solutions. These authors 


Rock em, 


264 Annals New York Academy of Sciences 


z= 4 Thermodynamic Model 


z= 4, Kinetic Model 
z=6, Kinetic Model 


z=8 
Kinetic Model 
z= 8, Thermodynamic 
Model 


0 0.2 0.4 0.6 0.8 1.0 
MOLAR FRACTION n 

FIGURE 7, Conditions for the zero net flux; comparison of the kinetic 
the thermodynamic models. 
derived a simple equation in differential form for the diffusion coetfidl: 
in concentrated solutions. They base their result upon the assumptil 
that it is the activity, and not the concentration, that is the driving for 
in diffusion. The expression for the diffusion coefficient is given 
these authors as follows: 


D d log a, ; ) 
ttt cast rl ats ail, 

Do dlogn 
where D is the diffusion coefficient of the concentrated solution, Do 
the diffusion coefficient of an ideal solution, and a, is the activity 
the species A. D/D, is what we calculated as the term proportional 


the net flux. 
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In order to compare our results with Equation III,1 we solve EQUA- 
pon III,5 for our thermodynamic model. This is done as follows: 
The activity of molecules A is defined as 


kT log a4 = pa — wo (11,6) 


here 4 — uA is the chemical potential of the solution with respect to 
} arbitrarily chosen standard solution whose chemical potential is As 


rom 
aap 
HA = ona ae 
e obtain 
d log a, UD Ng A : Ns 
{= =—* (WA +Ne) | —___ - 
d log n Do Ng Na (Na +Nz ) 
2Wap Ng 1 Ng 
(Na +N, ne 4N Gs: 
» Simply 
ze 
—_ = 1-2n(1-n) I, 
Do kT ities) 


We observe that the net flux is a quadratic function of n irrespective of 
while, in our kinetic treatment, D/D, is a polynomial in n of the order 
2(z- 1). 


= PART IV 
We assume now that we have 2 homogeneous binary solutions, of differ- 
concentrations, in contact. This means that n (x) is a discontinuous 
ion of x of the form shown in Figure 8. The concentration of mole- 
les of the species A to the left of the point of discontinuity x = x* is 
a to (n—d), while the concentration of the species A to the right of 
e point of discontinuity is equal to (a +d). 
Before proceeding with the derivation of the equation for the flux based 
FIGURE 8, let us pause and consider the various possibilities that 
‘model implies. 
The model allows each of the 2 phases that are in contact to be either 
a stable state, in a metastable state, or in an absolutely unstable 
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x* x— 


FIGURE 8, n (x) as a discontinuous function of x (see text), 


Let the molar free energy of phase 1 be F, and the molar free energy + 
phase 2 be F,. Let n, be the molar fraction of component A in phase 
and n, be the molar fraction of component A in phase 2. Then the con 
tion for the system to be in equilibrium with regard to the oxchanaan 
material across the boundary is 


dF, dF, | 

dn 1 7 dn 2 : 

At constant temperature and pressure, . 
dF, =(pt — py) dn, | 


and 


dF, = (u? —p2 ) dn, 


metastable state is 


Therefore, the condition for the binary solutions to form a stable or 
s 
ei —er =H2—w? | 
This is a necessary but not sufficient condition for an ‘‘equilibrium.’’ ¥ 
must still introduce limitations on each of the 2 solutions separately t 
prevent them from breaking spontaneously into 2 phases, even without : 
mutual contact. At constant temperature and pressure, we require that — 


d?F, d*F, 
—— >0 and >0 
dn,? dn,? 


If phases 1. and 2 are both stable, then these conditions are met if a 
only if 


ee 


Mazur & Alfrey: Supersaturated Binary Solutions 267 


‘hat is, the concentrations of component A will correspond to the 2 min- 
ma indicated by n, and n, (FiGuRE 1). However, if one of the 2 phases 
S in a metastable state, the entire system can form a metastable state 
ince, generally, pj“ + p/. 

Our task will be to investigate, with the aid of the kinetic model, what 
recisely are the conditions for the, zero flux and how they compare with 
esults of thermodynamic models, This is done again by employing a 
netic model for the derivation of the equation for the total net flux of 
omponent A across the boundary. In many details the derivation of an 
xpression for the total flux will be a repetition of our first problem. 

As before, we consider a molécular-pair interchange in:a-binary lattice 
FIGURE 9). However, the concentration of molecules to the left of the 
oundary is constant, and so -is the concentration of the molecules to the 
ight of the boundary (or point of discontinuity). The number of A-B inter 
hanges will be given by: (The number of sites per layer) - (the probabil- 
ty that a given pair of sites accommodates an A-B pair) - (the probability 
hat the atmosphere of neighbors is the particular one) « (e st/ eas 
eg bu ad ). A similar expression will’ be obtained for the namber of 
recisely reverse interchanges B-A, except that the second term denotes 
he probability that a given pair of sites accommodates a B-A pair, and 
he last term will be e"""?"*" , The net flux of A to the right. resulting 
rom this pair of interchanges is the difference between these 2 terms. | 


FIGURE 9, Molecular pair interchange in a binary lattice. 
: 
We shall now proceed to derive the individual terms in the expression 


or the net flux of A and then, by summing this expression over all pos- 
sible arrangements of molecules in the surroundings or atmospheres, we 
hall obtain an expression for the net flux of A to the right of the bound- 
ry point x = x*. 

The probability that the given pair of sites accommodates molecular 


: 
- 


a 
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pair AB is proportional to: 


D ao umn i laltematiear) dv,1 


The probability that the atmosphere of neighbors has a given arrangemé 
(described by N, and i, where N, is the number of molecules A in th 
atmosphere, while i is the number of molecules A in the atmos phere t 
the left of the point of the discontinuity x*) i is 


“Cc : at d)' (l-n+ a" at (n+ a)*4" nf ~n—d)*’ 2-NA # 

7 The net flux from the vite of interchianges-is praontiaiias to the ter 
C [Dap gE ru/kt arb, eo Bier /sT . ave 

where Dap *! given by EQUATION IV,1 and 
Dga = (l1—n+d) (n+d) ava 


Since, according to the adopted model for the energy of activation of the 
transition from state I to state II, 


‘ 
-E k = omA 
2 wor /kT A Pl /kT em e/2kT 


: 
-E kT = 
re rer/, ~e7% /kT o™ Ac/2xT 


we obtain the following expression for the net flux pees molecular inte ao 
change 


o=2K[ =! Q-a+gh/™ (n+ d) "48 (btn) 9 ee 


—d sh 


Let sinh BS, sen * CO mAs a 
| 2RT “2ET | 


The total flux will be bhathin by siimwttn EQUATION IV, 5 over 
possible arrangements of the molecules’ in® the Rewspnees with the fol 
lowing restrictions imposed on i and Ng’: 12e9 


= 


OSI NA; O<Na<N - 
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1e total flux is now 


N/2 
Total flux = 2K | (1—n) ee (’ ) @ >) 
: Pa : = i= 


avi 


Sere d_) (sty ( 24g?) a 
l-nid = n—-d vit N;0) 
ipa ie aa te Na- 2i 
sinh (M4 = 20 a i) Ae | —dcosh Wana) a2) Ae 
kT {s | KT 


nce m = 2(N, —2i). 


QUATION IV,6 is developed in APPENDIX II. Here, as in our first 
oblem, the total flux can be expressed in polynomial form as follows: 


(IV ,7) 


P 
_, pAe A 
petal flux = K > A, (n) sinh i — B, (n) cosh so 


aere A, (n) and B, (n) denote polynomials in n of order of 2 (z — 1). 


In spite of the apparent similarity of EQUATION IV,7 to EQUATION 
8 the polynomials involved are of entirely different form. The constant 
differs from the K in rquaTion II,8. However, as it was in the previ- 
is case for Ac > 0, the equation for the total flux is composed of 2 
ms, the first of which describes the contragradient part of the flow 
toss the boundary, while the second term describes the cogradient part 
the flow. When Ac <0, both Baye are oe” and the total flux is 
ways negative. 

Numerical computations of Equation A Il, 2 (APPENDIX I) were per- 
rmed on the electronic digital computer. For Ac /kT, 7 possible values 
are selected at successive intervals of, 0.2, covering the range of 
<Ac/kT < ‘1.4: Negative values of Ac/kT were not used, since then 
= flux will always be cogradient. When Ac/kT exceeds 1.5, the flux is 
ntragradient for almost every value of n, except for a very narrow strip 
ar the ‘‘edges’’ (n = 0 and n = 1). Successive values of both n and d 
sre increased by intervals An = 0. 05 and Ad = 0.05. Since the equation 
flux remains unchanged when n and 1—n are interchanged, only 
lues of n < 0.5 needed to be considered. 

The his of the computations are presented in FIGURES 10 to 13. 


Oa ny: 


€ 


FIGURE 10, The dependence of the net flux on the concentration ditterel 
at the a igtcephase boundary Ac! SE = 0.80. 
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FIGURE 11. The dependence of the net flux on the concentration diffe 
at the interphase boundary Ac /kT = 100. 
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FIGURE 12. The dependence of the net flux on the concentration difference 
| the interphase boundary Ae/kT = 1,20. 
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_ FIGURE 13, The dependence of the net flux on the concentration difference 
the interphase boundary Ae/kT = 1.40. - 
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For Ac/kT < 0.8 the net flux is always cogradient, and therefore we di 
not bother to plot these results. In rrcurEs 10 to 13, the total flux (on 
more correctly, the quantity proportional to the total flux) is plotted asi 
function of d for given values of n and Ac/kT. : 
The results of these computations throw light upon the problem of whe 
we call ‘‘metastable equilibria’? between 2 liquid phases of differen 
concentrations. In order to be able to interpret our results in terms o 
phase ‘‘pseudoequilibria,’’ it is important to examine the particular case¢ 
when the net flux is zero and to make use of a thermodynamic model. : 
TABLE 1, which is prepared from the data of rrcureEs 10, 11, 12, ana 
13, describes the conditions for zero flux. In this table, (nm + d) is tha 
concentration of phase 2 in thermodynamic equilibrium with the concens 
tration of phase 1, which is (nm — d). We notice that as (nm — d) is plotte: 


TABLE 1 


CONDITIONS FOR ZERO FLUX CONCENTRATION DIFFERENCE 
OF THE POINT OF DISCONTINUITY (=2D) 


Ae/kT =0.8 


versus (n + d) as in FicuReE 14, the curves have minima. These minin 

gradually move toward lower values of (n + d) and higher values of (n -d 
as Ac approaches Ac . By Ae we understand the smallest energy tha 
will still bring a contragradient flux. As Ac approaches A¢*, the point 0: 


have identical values of dF/dn. Suppose that we have phase 1, whe | 
composition in terms of molar fraction of component A can be depicted by 
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FIGURE 14. Conditions for zero net flux concentration of phase I versus 
icentration of phase II for several values of Ac /kT. 


point on an excess free-energy curve to the left of the point A (see 
GURE 15). This phase can be in an equilibrium with phase 2, of either 
2 possible compositions that are located to the right and to the left of 
+ inflection point DJ. One of them will correspond to the absolutely 
stable configuration (the part of the curve between points C and D), 
d the other to the metastable configuration (between points D and B). 
the concentration under consideration reaches point X, we obtain a 
se where the phase under consideration can form no such equilibrium 
h any other configuration to the right of point A. It is obvious th: t this 
nt X will correspond to the minimum value of (n — d) of Ficure 14 
{1 TaBLe 1. 

is situation can be described best on the hypothetical (dF/dn) 
sus n curve (FIGURE 16): the points X, A, B, C, and D are exactly 
| same points as indicated on ricureE 15. The part of the curve be- 
en points C and D represents absolutely unstable states. Again, we 
rt with point A. The corresponding point that represents phase 2 in an 
lilibrium with it is point B. As we move down to the left of point A, 
_ corresponding pseudoequilibrium concentration of phase B moves 
yn along the BD curve. As we reach point X, we reach the minimum 
ssible concentration of phase 1 that can be in contact with some con- 
tration corresponding to phase 2. 


so oh sieht 


0 0.5 1.0 
MOLAR FRACTION n 


FIGURE 15, Free energy of a binary mixture. The interconnected dots 


indicate pseudoequilibria between 2 binary mixtures. 


MOLAR FRACTION n 
FIGURE 16, Pseudoequilibria between 2 binary mixtures, 
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Going back to our kinetic model, we see that each concentration of 
tase 1 (indicated by [n —d]) above acertain value corresponds to 2 val- 
‘s of concentrations of phase 2 (indicated by [n + d]). 
Only one of these values is physically possible. These permissible 
mcentrations lie to the right of the minima of the (n — d) versus (n + d) 
irves (points between B and D of ricures 15 and 16). The concentra- 
ons to the left of these minima are physically impossible, since they 
ill fall into the absolutely unstable region between points C and D of 
e same figures (where either d?F/dn,? or d*F/dn,” <0). 
Thus far our kinetic model agrees qualitatively with the thermodynamic 
odel. In order to make a quantitative comparison, we calculate the free- 
lergy curve of binary solutions, assuming that the molecules are ran- 
ymly distributed (the so-called ‘‘zeroth-order’’ approximation of Bragg 
id Williams). The results are given in ricureE 17. In this figure 
F,/dn, is plotted against n,, for Ac/kT = 0.8. The calculations are 
used upon the equation 
Z£ sae log m4 + (1-2n,) 
1 pea Ely LF 


ze 


v,8) 


kT 


10) 0.5 1.0 
MOLAR FRACTION n 


FIGURE 17. Differences isi concentration (dF/dN) versus the concentration 
culated according to EQUATION V,10 for Ac /kT = 0.80. 
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In FIGURE 17, all points of interest (X, A, B, C, D) are indicated is 
identical fashion as in FiGuREs 15 and 16. From F1cuRE 17 we obtaiil 
the 2 concentrations (n — d) and (n + d),which are in zero-flux contaa 
(identical values of dF/dn); the results are plotted in FiGURE 14 (brokes 
line). We see that the shape of the (nm — d) versus (n + d) curves is very 
similar to that of the curves obtained from the kinetic model, except fo: 
the energy factor. The results derived from E QuaATION IV,8 for Ac/kT = 
0.8 are similar to the results based upon the kinetic model with Ac/kT 2 
1.15. In spite of the poor quantitative agreement, the thermodynamic ana 
kinetic treatments lead to the same qualitative conclusion, which may be 
stated as follows: : 

It is commonplace that, when 2 saturated solutions (m, and no ot 
FIGURE 1) are in contact, there is a zero diffusional flux across the 
interfacial boundary despite the discontinuity in concentration. It would 
now appear possible, in principle, to observe a zero flux across the 
boundary between 2 solutions that are not in equilibrium if the concentra 
tions are properly chosen. Compared with the saturated solutions (m, and 
n,), both of these phases must be overrich in the same species. Thus, 
one phase will be unsaturated, the other supersaturated. If the supersat- 
urated phase is metastable (rather than absolutely unstable), the 2-phase 
system as a whole would appear to be metastable. 


7 
' 
SUMMARY ; 
The diffusion behavior of supersaturated binary solutions is treate 
theoretically on the basis of a lattice model. A supersaturated solutior 
may be either metastable or absolutely unstable. In the region of absolute 
instability the diffusion coefficient is negative; fluctuations in concen 
tration are amplified by diffusion, and phase separation occurs spo. 
taneously. 
The theory also indicates the possibility of metastable 2-phase s 
tems in which 2 solutions in contact exhibit a zero diffusional flux across 


the interfacial surface even though the 2 concentrations are not in thermo- 


dynamic equilibrium. { 


APPENDIX I : 
In order to calculate the total flux of component A to the right oh 
positive x axis), we multiply Equations II,7 by the number of equiva 
lent configurations, and we then add these products for all possibl 
arrangements of the molecules in the atmosphere. This is done as follo 
Let the i molecules A be distributed among the sites at the left side o 
the atmosphere in such a way that jo of them will occupy the sites “" 


; 
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=0, and i — jy will occupy the sites with 7 = 1. Similarly, j2 of the 
maining N,4 — i molecules of kind A will occupy the sites with / = 3, 
id N4 —1—j, will occupy the sites with 1 = 2. Also, let a denote the 
umber of sites along 1 =0 and / =3, and N/2-— a denote the number of 
tes along / = 1 and / = 2. The total number of possible arrangements of 
4 molecules with given values of j, and j, and i among N sites is 


oie 2A 
2 (AI, 1) 
mae AM, -—1-js 


he energy of each possible arrangement depends always only upon 2 
ariables i and N,, since each molecule in the atmosphere forms the 
ime bond with the central molecule, irrespective of a. 

The equation for the total flux is now (from EquaTion II,7 and after 
troducing a new variable p = N, — 21): 


oa 
2 


Potal flux =k Ax © Ass De > > & (Na, i; jor 3) 


=-N 
PAO Poa. tA 9 j3=0 


. N. -p Naz +p ; ‘ 


Na -P Na +P rok Re 
cee aay 3 2. sinh - 
n 9 Jo ae a( 9 is) J3 kT S 
(AI,2) 


yme terms that obviously are odd functions of p have been dropped from 
euaTion AI,2. In Equation AI,2 


| ee Ne N\ [a 
6 N,, P; jor i3)= ‘ ; Na-P_., Nat+P _. 
Jo J3 ea 2 se 
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Therefore, 


N * a a 7 ! 
Total flux = K’ Pas re Pa = & (Na, P jo» is) 
e rae 


Na =0 at= fo=0.2 j3=0 


- {| 3N, +p he | 
7 ik ay) Be tay [sme z +i3~io}] 2 sinh ner hc/er ¢ 


’ Se ? 
K" = KAx — (Al,. 


Let us sum over positive and over negative values of p separately, ni 
such a way that the summation over p will run from 0 to W/2. Then 


N sone a ys } 
Total flux = K' og > ys oR nXA (1-1) *-"4@ 
, Na= P=1 = j9=0 4370 


2 
3N,4 - 
( m P sii) 2 sinh —— ( ~K' > > Na 
Na= p=t 
N-Na pAe 2 
(1 -n) & (Na, p) \ 2 cosh apron (Al,4) 


where parts I and II are enclosed by braces (j $). As we show below, the 
first part describes the contribution of the contragradient flow to the tota 
flux, while the second part describes the contribution of the cogradient 
flow toward the total flow when Ac >. This is so since part II is inde 
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dent of the arrangement of molecules inside the atmosphere. In 
UATION AI,4 


eee N a a 
: - > Sinko 
& (N,, P, Jo» is) = 2 
AES ABP < j j 
ree Alp TNS aee 
N = N 
=) — > —-a a a 
$(N.,-—D, rig l= 
Ng ep 5. Na 'P. os 


J3 
ou. N_ 
2 2 
g (V4 ? Pp) = (AI,5) 
Na+pP Na -P 
2 By 


now have an equation for the total flux in its final form. We will con- 
er the case of a= 1, only. After summing EquaTion AI,5 Gver jy and 
we obtain the following contributions of the contragradient and co- 
dient flows to the total flux: 


No ON ef | 
2 oF Sr 
K' Bit Slee birt ates * 
ge Na -P)\ Na +P 
i Na=0 p= 2 2 
N N art oS) 
ie aa oe 2 
+ (p+ - 
V4 —P Na +P Na-P _ Na +P 


Aes 
“iL aly § oe | 
2 2 pe 
— é/ - 2 sinh (AI, 6) 
(2 -p) kT 
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N = N N 
F a NA N-NA 2 2 
Il=K > > n “ (1—n) om tie 
y es A 
aye Dig Pty 


Na=0 p=i1 


N\ 2 
A ud Ne z j 
2 cosh ae + ae nA (lenin NA (AI,7 
Na 
Na =0 2. 


Note that when p = 0, J = 0 and JJ = 1, by binomial expansion (z QUATIONI 


AI, 7). 

Before summing parts I and II over p and N, we should stress the fac 
that the ‘factorial’? terms represented by symbols (*) are composed f 
integers only. Therefore, N, is an odd integer when p is an odd integ 
and N,4 is an even integer when p is an even integer. The final results 
for the total flux are: | 


CasEI: z=4 (N = 6): 
Total flux = K‘{n* (1 —n)* [10 sinh 3x — 2 cosh 3x] 
+[n* (1-—n)? +n? (1-n)*]- [20 sinh 2x -— 6 cosh 2x] 
~ [n> (l-n) + 3n3 (1-—n)? +n (1—n)5]- [10 sinh x — 6 cosh x] 
—[n® + 9n* (1-n)? + 9n? (1-n)* + (1-2)9]} (AI,: , 


CASE Il: z=6 (N = 10): 


Total flux = K‘ a® (1-n)5 [14 sinh 5x — 2 cosh 5x] 
+[n® (l—n)* +n* (1—7)®] [56 sinh 4x — 10 cosh 4x] 
+ [n7 (1 -n)3 + 2.5n5 (1 —n)5 + n3 (1 —2n)7] [84 sinh 3x — 20 cosh 3x | 
+ [n® (l-n)? + 5n® (l—n)4 +5n* (1—n)® +n? (1—n)8] 
[56 sinh 2x — 20 cosh 2x] + [n9 (l—n) + 10n’7 (1—n)? 
+ 20n° (1 -n)5+ 10n3 (1=n)74n(1 ~n)°] [14 sinh x -10 cosh x] 
—[n® + 25n® (1 —n)? + 100n® (1 ~n)* + 100n4 (1 —n)® 


+25n? (1—n)®§ +(1-—n)” (AL, 
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ASE III: z=8 (NW = 14): 


Total flux=K' n’7 (1—n)’7 [18 sinh 7x — 2 cosh 7x] 
+ [n® (1 —n)® + n° (1 —n)®] [108 sinh 6x — 14 cosh 6x] 
fr? (1 n)5 + Sn7(1— )? 405 (1 = n)°] [270 sinh § x—42 cosh 5 x] 
+ et =n) + 4.2n*® (1 —n)° 4 4.2n® (1 —n)® +'n® (1-2) ”] 
[360 sinh 4 x — 70 cosh 4 x] + [n"1(1 — n)* + 7n? (1 —n)5 
12.6n’7 (1—n)’ + 7n5(1—n)? + n3(1—n) "] [270 sinh 3x — 70 cosh 3 x] 
+ [a (1 —n)? + > 2 (1 ~ n)* + 358 (1 ~n)® + 3509 (1 —a)® 
+2 n* (1~n)” +n? (1-n)”]- [108 sinh 2x — 42 cosh 2x] 
4 [a (1-1) + 21n™ (1-2)? 10509 (1 0)5 + 17507 (1-2)? 
asa -0)°4: Waeil ~n)2 He(tan)) 
[18 sinh x — 14 cosh x] —[n™ + 49n (1—n)? + 441n”™ (1 -n)* 
+ 1225n® (1 —n)® 1225n® (1 —n)® + 441n* (1—n)” 
+ 49n? (l—n)?+(1—n)*] (AI,10) 
‘all these equations x =Ac¢/kT. 
| APPENDIX II 


We rewrite & QUATION IV,6 by introducing a new variable p = N, — 2i. 


Total flux = 2K | (1—n)? — d? 


dg as 
= 
Fas 
ae Br 
N+ 
ae) 
=) ~ 
17s 
ae | 
+1) Q& 
A. 
> 4 
sf 
a] 
food ‘\ 
) is} 
ml + 
1] a 
Q. 
ee 
SES 
v 


pe pe 
— n2 +d?) sinh — dcosh (All, 1) 
(n — n? +d?) sin = a | 


hb debe 244 we i 
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The restrictions on p and Ng are: 


-N ¥ N 
we =e Pees 
eer ae 
Here again, both p and Ng can be only odd integers or even inte 
simultaneously. In order to make EQUATION 1 more compact, the fo 
ing abbreviations are introduced: 


n+d n-d 


eae eae Si SS, 


Then for the specific case of N = 6 we obtain 


Total flux = 2K [(1 —n)?- d?}° {c — n* + d?) G = B?) sinh = 


Ae 


Dy) « 
+ 3 (a?~ 8?) (1 + a) sinh —— + 3 (a— f) (1 + 308 + a7p”) sh 


r 
« 


3Ae 
-—d |(a* + B%) cosh +3 (a? + B?) (1 q 
[ a“ +B + af) (AlL.2 
2Ae Ae 
cosh + 3(a0 B) (1+ 3a + a?6*) cosh — 


kT 
+ (1+ 9aB + 9a7B? + <p} 


For n = 0.5, aB = 1, and we obtain 
3Ae 


Total flux = 2K (0.5 — d)® {02s +d?) G — B*) sinh a 


2Ae Ae 
+ 6 (a? - B?) sinh —— 15 (a - baits coreg 
B*) sin + (a — B) sinh = 


Ae 


2A 
+ 6 (a? + B?) cosh 


-~d G + B%) cosh S 


Ae 
+ 14 (a0 B) cosh a + 20) } (AIL, 
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In the limiting case of Ae = 0 we have 

Total flux = -2Kd (AII,4) 
ace (1 — n)* — d? = (10 a) (10 8) and from binomial expansion of 
QuaTion AII,2. For the case in which pure solutions are in contact, 


=d, and we have 


eArar loeb 
< Total flux=-2Kn |1-2n(l-e ) <0] PAIS) 


” 


Forn=0.5 andn=d we have 


Total flux = —Ke “4/*7 (All,6) 
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